Absence of epidemic threshold in scale-free networks with connectivity correlations 
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Random scale-free networks have the pecuhar property of being prone to the spreading of in- 
fections. Here we provide an exact result showing that a scale-free connectivity distribution with 
diverging second moment is a sufficient condition to have null epidemic threshold in unstructured 
networks with either assortative or dissortative mixing. Connectivity correlations result therefore 
ininfluential for the epidemic spreading picture in these scale-free networks. The present result is 
related to the divergence of the average nearest neighbors connectivity, enforced by the connectivity 
detailed balance condition. 
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Complex networks play a capital role in the model- 
ing of many social, natural, and technological systems 
which are characterized by peculiar topological proper- 
ties [0, ||. In particular, small- world properties ||] and 
scale-free connectivity distributions Q appear as com- 
mon features of many real- world networks. The statis- 
tical physics approach has been proved a very valuable 
tool for the study of these networks, and several surpris- 
ing results concerning dynamical processes taking place 
on complex networks have been repeatedly reported. In 
particular, the absence of the percolation l^j, ^ and epi- 
demic [0, |[ ^, |l^ thresholds in scale-free (SF) networks 
has hit the community because of its potential practical 
implications. The absence of the percolative threshold, 
indeed, prompts to an exceptional tolerance to random 
damages On the other hand, the lack of any epi- 

demic threshold makes SF networks the ideal media for 
the propagation of infections, bugs, or unsolicited infor- 
mations 01 . 

Recent studies have focused in a more detailed topo- 
logical characterization of several social and technological 
networks. In particular, it has been recognized that many 
of these networks possess, along with SF properties, non- 
trivial connectivity correlations |p^ . For instance, many 
social networks show that nodes with high connectivity 
will connect more preferably to highly connected nodes 
; a property referred to as "assortative mixing" . On 
the opposite side, many technological and biological net- 
works show "dissortative mixing" ; i. e. highly connected 
nodes are preferably connected to nodes with low con- 
nectivity ^ Correlations are very important 
in determining the physical properties of these networks 
and several recent works are addressing the effect of 
"dissortative mixing" correlations in epidemic spreading 
| p7| , p^ , p^ . The fact that highly connected nodes (hubs) 
are more likely to transmit the infection to poorly con- 
nected nodes could somehow slow down the spreading 
process. By numerical simulations and analytical argu- 
ments it has been claimed that, if strong enough, connec- 



tivity correlations might reintroduce an epidemic thresh- 
old in SF networks, thus restoring the standard tolerance 
to infections. 

In this paper we analyze in detail the conditions for the 
lack of an epidemic threshold in the susceptible-infected- 
susceptiblc model in SF networks. We find the exact 
result that a SF connectivity distribution P{k) ~ k^"' 
with 2 < 7 < 3 m unstructured networks with assorta- 
tive or dissortative mixing is a sufficient condition for a 
null epidemic threshold in the thermodynamic limit. In 
other words, the presence of two-point connectivity cor- 
relations does not alter the extreme weakness of SF net- 
works to epidemic diffusion. This result is related to the 
divergence of the nearest neighbors average connectivity, 
divergence that is ensured by the connectivity detailed 
balance condition, to be satisfied in physical networks. 
The present analysis can be easily generalized to more 
sophisticated epidemic models. 

In the following we shall consider unstructured undi- 
rected SF networks, in which the only relevant property 
of the nodes is their connectivity [ pT| , with distribution 
P{k) - Ck-^, with 2 < 7 < 3. P(fc) is defined as the 
probability that a randomly selected node has k connec- 
tions to other nodes. In this case the network has un- 
bounded connectivity fluctuations, signalled by a diverg- 
ing second moment (k!^) — > oo in the thermodynamic 
limit fee — > oo, where k^ is the maximum connectivity of 
the network. It is worth recalling that in growing net- 
works kc is related to the network size N askc '-^ ]\ri/(7-i) 
1^. Finally, we shall consider that the network presents 
assortative or dissortative mixing allowing for non triv- 
ial two-point connectivity correlations. This corresponds 
to allow a general form for the conditional probability, 
P{k' I fc), that a link emanated by a node of connectivity 
k points to a node of connectivity k' . 

As a prototypical example for examining the proper- 
ties of epidemic dynamics in SF networks we consider 
the susceptible-infected-susceptible (SIS) model in 
which each node represents an individual of the popu- 
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lation and the links represent the physical interactions 
among which the infection propagates. Each individual 
can be either in a susceptible or infected state. Suscep- 
tible individuals become infected with probability A if at 
least one of the neighbors is infected. Infected nodes, 
on the other hand, recover and become susceptible again 
with probability one. A different recovery probability can 
be considered by a proper rescaling of A and the time. 
This model is conceived for representing endemic infec- 
tions which do not confer permanent immunity, allowing 
individuals to go through the stochastic cycle suscepti- 
ble — > infected — > susceptible by contracting the infection 
over and over again. In regular homogeneous networks, 
in which each node has more or less the same number of 
connections, fc ~ {k), it is possible to understand the be- 
havior of the model by looking at the average density of 
infected individuals p{t) (the prevalence). It is found that 
for a spreading probability A > Ac, where Ac is the epi- 
demic threshold depending on the network average con- 
nectivity and topology, the system reaches an endemic 
state with a finite stationary density p. If A < Ac, the 
system falls in a finite time in a healthy state with no 
infected individuals (p = 0) . 

In SF networks the average connectivity is highly fluc- 
tuating and the approximation k ~ (k) is totally inade- 
quate. To take into account the effect of the connectivity 
fluctuations, it has been shown that it is appropriate to 
consider the quantity pk |0, ||, ^ , defined as the density 
of infected nodes within each connectivity class k. This 
description assumes that the network is unstructured and 
that the classification of nodes according only to their 
connectivity is meaningful |21 . Following Ref. iQ, the 
mean-field rate equations describing the system can be 
written as 



dpkjt) 
dt 



-Pk{t) + Afc [1 - pk{t)] J2 P{k' I k)pu\t). (1) 



The first term on the r.h.s. represents the annihilation 
of infected individuals due to recovery with unitary rate. 
The creation term is proportional to the density of sus- 
ceptible individuals, 1 — pk, times the spreading rate, 
A, the number of neighboring nodes, k, and the prob- 
ability that any neighboring node is infected. The lat- 
ter is the average over all connectivities of the prob- 
ability P{k' I k)pki that a link emanated from a node 
with connectivity k points to an infected node with con- 
nectivity k' . It is worth remarking that, while keeping 
into account the two point connectivity correlations, as 
given by the conditional probability P{k' \ k), yet we have 
neglected higher order density-density and connectivity 
correlations. Eq. (|l]) is therefore exact for the class of 
Markovian networks p^ , in the limit of low prevalence 
{pit) « 1). 

In the case of uncorrelated networks each link points, 
with probability proportional to k'P{k'), to a node of 



connectivity k' , regardless of the emanating node's con- 
nectivity. In this case, in the stationary state {dtp = 
0), J2k' ^ik' \k)pk'{t) assumes a constant value inde- 
pendent on k and t and the system can be solved 
self-consistently obtaining that the epidemic threshold is 
given by pi 



Ac 



{k^Y 



(2) 



For infinite SF networks with 7 < 3, we have Ik"^) — 
00, and correspondingly Ac = 0; i.e. uncorrelated SF 
networks allow a finite prevalence whatever the spreading 
rate A of the infection. Finally, from the solution of pk, 
one can compute the total prevalence p using the relation 

In the case of correlated networks the explicit solution 
of Eq. (0) is not generally accessible. However, it has 



been shown that the epidemic threshold is given by |16 



Ac 



A„ 



(3) 



where A^ is the largest eigenvalue of the connectivity 
matrix C, defined by Ckk' = kP{k'\k). In Ref. ^ it 
has been shown how this general formalism recovers pre- 
vious results for uncorrelated networks, obtaining that, 
in this case, A^ = {k'^) / (fc). More generally, by looking 
at Eq. (H), the absence of an epidemic threshold corre- 
sponds to a divergence of the largest eigenvalue of the 
connectivity matrix C in the limit of an infinite network 
size N —>■ 00. In order to provide some general state- 
ment on the conditions for such a divergence we can make 
use of the Frobenius theorem for non-negative irreducible 
matrices |22 . This theorem states the existence of the 
largest eigenvalue of any non-negative irreducible matrix, 
eigenvalue which is simple, positive, and has a positive 
eigenvector. One of the consequences of the theorem is 
that it provides a bound to such largest eigenvalue [ p3[ . 
In our case the matrix of interest is the connectivity ma- 
trix and, since C is non- negative and irreducible |24[| , it 
is possible to find lower and upper bounds of A™. In 
particular, we can write ||23| 



A^>min^5]fc'£P(£|fc)P(fc'K). 



(4) 



This inequality relates the lower bound of the largest 
eigenvalue A^ to the connectivity correlation function 
and, as we shall see, allows to find a sufficient condition 
for the absence of the epidemic threshold. 

In order to provide an explicit bound to the largest 
eigenvalue we must exploit the properties of the condi- 
tional probability P{k' | fc). A key relation holding for all 
physical networks is that all links must point from one 
node to another. This is translated in the connectivity 
detailed balance condition iflql 



kP{k' I k)P{k) = k'P{k I k')P{k') 



(5) 
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which states that the total number of Hnks pointing from 
nodes with connectivity k to nodes of connectivity k' 
must be equal to the total number of links that point 
from nodes with connectivity k' to nodes of connectiv- 
ity k. This relation is extremely important since it con- 
straints the possible form of the conditional probability 
P{k' I k) once P{k) is given. By multiplying by a fc fac- 
tor both terms of Eq. and summing over fc' and fc, we 
obtain 

(fc2)^^fc'P(fc')^A:P(fc|fc'), (6) 

fe' k 

where we have used the normalization conditions 
Efe^(fc) = Efc'^(^'l^) = 1- The term fc„„(fc',fcc) = 
fcP(fc I fc') defines the average nearest neighbor con- 
nectivity (ANNC) of nodes of connectivity fc'. This is 
a quantity customarily measured in SF and complex 
networks in order to quantify degree-degree correlations 
0, |l|, The dependence on fcc is originated by the 
upper cut-off of the fc-sum and it must be taken into ac- 
count since it is a possible source of divergences in the 
thermodynamic limit. In SF networks with 2 < 7 < 3 
we have that the second moment of the connectivity dis- 
tribution diverges as (fc^) k^~^ We thus obtain 
that 

J2 k'P{k')Kn{k\ fcc) ^ TTT^fcc"^- (7) 

^ (3-7) 

In the case of dissortative mixing ||l^, the function 
fcnn(fc',fcc) is decreasing with fc' and, since fc'P(fc') is an 
integrable function, the l.h.s. of Eq. (Q) has no diver- 
gence related to the sum over fc'. This implies that the 
divergence must be contained in the fcc dependence of 
knn{k' , kc). In other words, the function knn{k' , kc) — > 00 
for fcc 00 in a non-zero measure set. In the case of as- 
sortative mixing, fcnn(fc ; fcc) is an increasing function of 
fc' and, depending on its rate of growth, there may be sin- 
gularities associated to the sum over fc'. Therefore, this 
case has to be analyzed in detail. Let us assume that 
the ANNC grows as fcnn(fc',fcc) ^ ak'f^, (3 > 0, when 
fc'— i'CxD. If/3<7 — 2, again there is no singularity re- 
lated to the sum over fc' and the previous argument for 
dissortative mixing holds. When 7 — 2 < /3 < 1 there 
is a singularity coming from the sum over fc' of the type 
afcf~^^~^\ However, since Eq. (Q) comes from an iden- 
tity, the singularity on the l.h.s. must match both the 
exponent of fcc and the prefactor on the r.h.s. In the case 
7 — 2 < /3 < 1, the singularity coming from the sum is 
not strong enough to match the r.h.s. of Eq. (0) since 
/3 — (7 — 2) < 3 — 7. Thus, the function fc„„(fc', fcc) must 
also diverge when fcc ^ 00 in a non-zero measure set. Fi- 
nally, when f3 > 1 the singularity associated to the sum 
is too strong, forcing the prefactor to scale as a ~ rk],~^ 
and the ANNC as fc„„(fc',fcc) ~ rkl~^k'l^. It is easy to 
realize that r < 1 , since the ANNC cannot be larger than 



fcc. Plugging the fc„„(fc', fcc) dependence into Eq. (|^) and 
simplifying common factors, we obtain the identity at the 
level of prefactors 



2-7 + /3 3-7' ^ ' 

Since (3 > 1 and r < 1, the prefactor in the l.h.s. of 
Eq. (Q) is smaller than the one of the r.h.s. This fact 
implies that the tail of the distribution in the l.h.s. of 
Eq. (|^) cannot account for the whole divergence of its 
r.h.s. This means that the sum is not the only source of 
divergences and, therefore, the ANNC must diverge at 
some other point pGf . 

The large fcc behavior of the ANNC can be plugged in 
Eq. (^) obtaining that 

A,2„ > mm ^ IP{1 1 fc)fc„„(^, fcc) (9) 

The r.h.s. of this equation is a sum of positive terms and 
diverges with fcc at least as fc„n(f , fcc) both in the dissor- 
tative or assortative cases . This readily implies that 
A„i > 00 for all networks with diverging (fc2). Finally 
Eq. (H) yields that the epidemic threshold vanishes in the 
thermodynamic limit in all SF networks with assortative 
and dissortative mixing if the connectivity distribution 
has a diverging second moment; i.e. a SF connectiv- 
ity distribution with exponent 2 < 7 < 3 is a sufficient 
condition for the absence of an epidemic threshold in un- 
structured networks with arbitrary two-point connectivity 
correlation function. 

In physical terms, the absence of the epidemic thresh- 
old is related to the divergence of the average nearest 
neighbors connectivity (knn) in SF networks. This 
function is defined by 

(fcnn)^ = P{k)knn{k, kc), (10) 

fc 

where we have explicitly considered fcc as a growing func- 
tion of the network size N. By using the analysis shown 
previously it follows that {knn) — * oo when iV — > oo. 
In SF networks this parameter takes into account the 
level of connectivity fluctuations and appears as ruling 
the epidemic spreading dynamics. Somehow the number 
of neighbors that can be infected in successive steps is 
the relevant quantity. Only in homogeneous networks, 
where (knn) ^ — (fc), the epidemic spreading properties 
can be related to the average connectivity. Noticeably, 
the power-law behavior of SF networks imposes a diver- 
gence of (fcrin)^ independently of the level of correlations 
present in the network. This amounts to lower to zero 
the epidemic threshold. On the practical side, connec- 
tivity correlation functions can be measured in several 
networks and show assortative or dissortative behavior 
depending on the system. These measurements are al- 
ways performed in the presence of a finite fcc that allows 
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the regularization of the function knn{k,kc)- The most 
convenient way to exploit the infinite size singularity is 
to measure the average nearest neighbor connectivity for 
incr paeing network eigoc. All SF notworkc with 2 < j < 3 
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ing in correlated complex networks (2002), e-print cond- 
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in structured scale-free networks (2002), e-print cond 
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ensured by the connectivity detailed balance condition 
alone. Thus it is a very reneral results boldinp; for all 



C.-H. Chang, and T. Krueger, Epidemic 
on scale-free graphs: the interplay b etween 



SF Networks with 2 < •y < ?,. On the contrary, the SF 
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behavior with 2 < 7 < 3 is a necessary condition for the 
lack of epidemic threshold only in networks with gen- 
eral two-point connectivity correlations and in absence 
of higher-order correlations. The reason is that the rela- 
tion between the epidemic threshold and the maximum 
eigenvalue of the connectivity matrix only holds for these 
classes of networks. Higher order correlations, or the 
presence of an underlying metric in the network [pTf , can 
modify the rate equation at the basis of the SIS model 
and may invalidate the present discussion. 
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From the Frobenius theorem [^2| it can be proved that 
the maximum eigenvalue. Am, of any non-negative ir- 
reducible matrix, A^k' , satisfies the inequality Am > 



where ip{k) is any positive vec- 
and «/'(fc) = fc we recover the 



tor. By setting A = 
inequality of Eq. (^ . 
The irreducible property of the connectivity matrix is a 
simple consequence of the fact that all the connectivity 
classes in the network are accessible. That is, starting 
from the connectivity class k it is always possible to find 
a path of links that connects this class to any other class 
k' of the network. If this is not the case it means that the 
network is built up of disconnected irreducible subnet- 
works and, therefore, we can apply the same line of rea- 
soning to each subnetwork. Notice that being irreducible 
is not equivalent to being fully connected at the node to 
node level, but at the class to class level. 
For 7 = 3 the second moment diverges as 
but the argument, though more involved, is still valid. 
The case (3 = 1, that is, k„n{k' , kc) ~ ak' when k' — > 00, 
is rather pathological because is the only case where the 
ANNC can be a non diverging function. If a < 1, Eq. (^) 
implies that it must exist at least one diverging point 
out of the one coming from the sum over k' and the dis- 
cussion made in the text is valid. However, when a = 1 
this divergence is enough to fulfill the identity (Q) and 
knn{k,kc) can be a convergent function when fee — > 00. 
Even in this pathological situation it is possible to prove 
that the maximum eigenvalue diverges in the thermody- 
namic limit. In order to do so the argumentation is to 
be translated to the second moment of the nearest neigh- 
bor's connectivity. Since the inequality {x'^) > (x)^ is 
true for any random variable, the function fc2„„(A;, kc) is 
either ~ k'^ , in which case the maximum eigenvalue is 
trivially 00, or the probability P{k' \ k) approaches to SF 
distribution when kc ^ 00 and Am diverges as well. 
One may argue that, since we are calculating a minimum 
for k, if the transition probability P{£\ko) is zero at some 
point ko, this minimum is zero. In this case it is possi- 
ble to show that repeating the same argument with C'^ 
instead of provides us an inequality that avoids this 
problem. 



